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Key words magnetohydrodynamics (MHD) - magnetic fields - turbulence 

The components of the total stress tensor (Reynolds stress plus Maxwell stress) are computed within the quasilinear 
approximation for a driven turbulence influenced by a large-scale magnetic background field. The conducting fluid has 
an arbitrary magnetic Prandtl number and the turbulence without the background field is assumed as homogeneous and 
isotropic with a free Strouhal number St. The total large-scale magnetic tension is always reduced by the turbulence with 
the possibility of a 'catastrophic quenching' for large magnetic Reynolds number Rm so that even its sign is reversed. The 
total magnetic pressure is enhanced by turbulence in the high-conductivity limit but it is reduced in the low-conductivity 
limit. Also in this case the sign of the total pressure may reverse but only for special turbulences with sufficiently large 
St > 1. The turbulence-induced terms of the stress tensor are suppressed by strong magnetic fields. For the tension term 
this quenching grows with the square of the Hartmann number of the magnetic field. For microscopic (i.e. small) diffu- 
sivity values the magnetic tension term becomes thus highly quenched even for field amplitudes much smaller than their 
equipartition value. In the opposite case of large-eddy simulations the magnetic quenching is only mild but then also the 
turbulence-induced Maxwell tensor components for weak fields remain rather small. 
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1 Introduction 

Differential rotation and fossil fields do not coexist. A 
nonuniform rotation law induces azimuthal fields SB^ from 
an original poloidal field Br which together transport an- 
gular momentum in radial direction reducing the shear Sfl 
via the large-scale Lorentz force J x £?, i.e. 



R 



sn _ BrSb^ 

St fiopAR 



(1) 



As the induced B^ results as 8B$ ~ AflBj^St the duration 
of the complete decay of the shear (i.e. <5fi = Ail) is St ~ 
Jpi^pR/ Br. This is a short time of order 10000 yr for a 
fossil field of 1 Gauss compared with the time scale of the 
star formation. All protostars should thus rotate rigidly. 

Equation (Q3 is also used for the explanation of the ob- 
served torsional oscillations of the Sun. With Br ~ 5 Gauss 
and Bcf, ~ 10 000 Gauss the estimation for RSfl is 10 m/s 
which is close to the observed value of 5 ms _1 . The result 
is that - if Eq. (Q]i is correct - the maximal field strength of 
the invisible toroidal fields should not be much higher than 
10000 Gauss. However, the solar convection zone is turbu- 
lent and it is not yet clear whether Eq. ([]]) is also true for 
conducting fluids with fluctuating flows and fields. 

In this paper the total Maxwell stress is thus derived 
for a turbulent fluid under the presence of a uniform back- 
ground field B. The fluctuating flow components are de- 



noted by u and the fluctuating field components are denoted 
by b. The standard Maxwell tensor 



Mi, 



1 



— BjBi 



-B 2 5, 



(2) 



1 

Mo " 2^o 
for the considered MHD turbulence turns into the general 
ized stress tensor 

T 



M}° = M,, 



pQ 



with the one-point correlation tensor 

Qij = (u t (x,t)uj(x,t)) 

of the flow and the turbulence-induced Maxwell tensor 



(3) 



(4) 



1 



1 



Mfj = —(h{ Xl t)b 3 {x,t)) - —(b'(x,t))6, 
Mo ^Mo 

The generalized Lorentz force F is then 



F t = M„ 



13,3 ' 



(5) 



(6) 



If the only preferred direction in the turbulence is the uni- 
form background field B both the tensors Qij and have 
the same form as the Maxwell tensor (|2|i but with two un- 
known scalar parameters. It makes thus sense to write 



M tot = _ K)B . B . _ J_ (1 _ Kp )B 2 S l3 



(7) 
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M0 ^M0 

for the total stress tensor ©. The first term of the RHS de- 
scribes a tension along the magnetic field lines while the 
second term is the sum of the magnetic-induced pressures 
transverse to the lines of force. The main role of the first 
term in stellar physics is an outward-directed angular mo- 
mentum transport if Br B^ < 0. If its coefficient 1 — k 
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would change its sign under the presence of turbulence then 
for the same magnetic geometry the angular momentum 
transport would be inwardly directed. The Lorentz force © 
with <[Vj> becomes 



(8) 



F=(1-k)J x B - ^—(k - k p ) VB 2 , 

so that the 'laminar' Lorentz force J x B has to be multi- 
plied with the factor 1 — k and an extra magnetic pressure 
appears if the n's are unequal (and they are) due to the ac- 
tion of the turbulence. If the k is positive then its ampli- 
tude should not exceed unity as otherwise the direction of 
the Lorentz force reversed. Roberts & Soward (1975) con- 
sidering only the terms of the Maxwell stress found (large) 
positive k and negative k p , i.e. k — — k p = rfr/v with t]t 
the well-known eddy diffusivity (see Eq. d42i i. below). 

Riidiger et al. (1986) found k also as positive and as 
running with the magnetic Reynolds number Rm of the tur- 
bulence even for Rm > 1. Kleeorin et al. (1989) suggest 
that K p > and even larger than unity so that the total mag- 
netic pressure changes its sign and becomes negative. The 
resulting instability may produce structures of concentrated 
magnetic field and may be important for sunspot formation 
(Kleeorin et al. 1990; Brandenburg et al. 2010, 201 1). 

Hence, the k's have an important physical meaning. In 
the simplest case they both would result as negative. Then 
the effective pressure is increased by the magnetic terms 
and also the tension term 1 — k remains positive so that 
the Lorentz force in turbulent media is simply amplified. 
Many more serious consequences would result from posi- 
tive k and K p if exceeding unity. In this case the Lorentz 
force changes its sign with dramatic consequences for the 
theory of torsional oscillations and solar oscillations (Klee- 
orin & Rogachevskii 1994: Kleeorin et al. 1996). 

The turbulence-induced modification of the Maxwell 
stress can also be important in other constellations where 
large-scale fields and turbulence simultaneously exist. 
Tachocline theory, jet theory, the structure of magnetized 
galactic disks (Battaner & Florido 1995) or oscillations of 
convective stars with magnetic fields could be mentioned. 



2 Equations 

We apply the quasilinear approximation known also as the 
second order correlation approximation (SOCA). All pre- 
liminary steps to derive main equations for the problem at 
hand were described by Riidiger & Kitchatinov (1990). The 
fluctuating magnetic and velocity fields are related by the 
equation 

b(k,w)= l{k ' B) u(k,u), (9) 
— iui + rjk z 

where the hat notation marks Fourier amplitudes, e.g. 
b(r,t)= { e i{kx ~ ut) b{k,Lu) dfcdw, (10) 



and the same for the velocity field. The influence of mean 
magnetic field on turbulence is described by the relation 



ii(k, u) = 



(fc,w) 



i | (k;Vr 



(ID 



where V = B/^/fiop is the Alfven velocity of the large- 
scale background field, rj is the microscopic magnetic dif- 
fusivity, and v is the microscopic viscosity. In Eq. ( fTTT ) ii 
is the (Fourier-transformed) velocity field modified by the 
mean magnetic field and stands for the velocity of the 
'original' turbulence which is assumed to exist for B = 0. 
The original turbulence is assumed as statistically homoge- 
neous and isotropic, i.e. 

^\k^)uf\k'M)) = E{k ^ 



16Trk 2 V k 2 
x 5(k + k')S(uj+uj'), (12) 

where E(k,uj) is the positive-definite spectrum function of 
the turbulence. Here 



u 2 = 



J J E(k,uj) dfcdw 



(13) 



defines the rms velocity u of the original turbulence. Equa- 
tions (0 to ( fT2l suffice to derive the values of the k's. 

3 Weak field 

We proceed by considering special cases. For weak mean 
magnetic field one finds from the expressions given by 
Riidiger & Kitchatinov (1990) the relation 



1 

15 



£fc 2 0(2?7 + j/)/c 4 -co 2 ) 
{v 2 k A +Lo 2 ){r 1 2 k i + oj 2 ) 



d/cdw, 



1 

15 



o o 



Ek 2 (v(8r]-is)k 4 -9w 2 ) 
(^ 2 fc 4 + u; 2 )(?7 2 fc 4 +w 2 ) 



dkduj. 



(14) 



The expressions do not have definite signs so that it remains 
unclear whether the large-scale Maxwell stress is increased 
or decreased by the turbulence. Even the signs of k and k p 
may depend on the spectrum of the turbulence. 

The simplest case is a turbulence with a white-noise 
spectrum containing all frequencies with the same ampli- 
tude. Here and in the following we shall use the Strouhal 
number St and the normalized characteristic frequency w* 



St = — T c , 



n 



(15) 



(l c correlation length, r c correlation time). The turbulence 
frequency w* measures the characteristic frequency of the 
turbulence spectrum in relation to the diffusion frequency. It 
is large for flat spectra such as white noise and it is small for 
very steep spectra like S functions. On the other hand, it is 
large in the high-conductivity limit and it is small in the low- 
conductivity limit. E.g., it is much larger than unity if the 
microscopic (Spitzer) diffusivity is used (high-conductivity 
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limit). It should be unity if - as it is used in large-eddy sim- 
ulations - rj = r)T — wl%. In the numerical integrations 
presented below the limit w* — > (i.e. low-conductivity 
limit) applies to the case of the frequency spectrum as a 
Dirac delta function S(ui). 

The product of St and w* giving the magnetic Reynolds 
number 

ul c 

Rm = — , (16) 

V 

where we have used the relation r c = 1/w as a definition of 
the correlation time. Then it is w* = Rm/St. 

In the high-conductivity limit ('white noise') one finds 
the simple results 

K = — Rm St, n v — Rm St, (17) 

15 p 15 

so that the k is positive and runs with Rm St = u 2 T c /rj 
which is the (large) ratio of the eddy diffusivity and the mi- 
croscopic diffusivity hence the magnetic tension is always 
(strongly) reduced. On the other hand, the magnetic pres- 
sure is increased (see Eq. [7J. The negative sign of the value 
of k p excludes the possibility that the effective pressure 
term in Eq. (0 changes its sign so that the total magnetic 
pressure becomes negative. This is formally possible after 
dTTb for the magnetic tension parameter 1 — k. 

The white-noise approximation, however, is not perfect. 
If, for example, the opposite frequency profile for very long 
correlation times, i.e. E oc 8{uS), is applied to (TBI then 
again the k is positive but the sign of n p depends on the 
magnetic Prandtl number 

Pm=-. (18) 
V 

It is thus necessary to discuss the integrals in (fT~4T > in more 
detail. 

3.1 Pm > 1 

For Pm > 8 one finds that n p is negative-definite for all 
possible spectral functions. The coefficient 1 — k p of the 
magnetic pressure is thus positive-definite and cannot be- 
come negative. This is not true for k. We shall show that the 
k will 'almost always' be positive so that the Lorentz force 
term in the generalized Lorentz force expression ([8]l is 'al- 
most always' quenched by the existence of the turbulence. 
For v = r) the expressions ( TT4b turn into 



1 

15 



Ek 2 {3f] 2 k 



2 7,4 



U! 



(u 2 +rj 2 k 4 ) 2 



oooo 



dk dui, 



Ek 2 (7r) 2 k 4 -9a; 2 

15 J J (w 2 + ?7 2 fc 4 ) 2 
o o 



dk dui, 



(19) 



which again do not have definite signs. One can write, how- 
ever, the expression for n as 



1 

15 



2?fk b E 



o o 



(w 2 + 7? 2 fc 4 ) 



-dkduj 



1 

15 



dE 



o o 



lo 2 + rfk 4 duj 



dfcdc 



(20) 



from which k proves to be positive-definite for all spectral 
functions E which do not increase for increasing to. We 
shall see that the positivity of k which reduces the effec- 
tivity of the angular momentum transport is a general result 
of the SOCA theory. 

Further simplifications can be achieved by applying the 
model spectrum 

2w 

E(k,u)=q{k) — ; ^ (21) 



7r(u> 2 + CJ 2 ) 



with 



oc 

J q(k) dk = u 2 , 



(22) 



where w is a characteristic frequency of the turbulence spec- 
trum. For w — > (f2Tb represents a Dirac <5-function while 
w — > oo gives 'white noise'. The results are 



w + Zrjk 2 



1 

15ry J (w + rjk 2 ) 2 
o 



and 



7rjk 2 



15?y J (w + rjk 2 ) 2 
o 



qdk 



qdk. 



(23) 



(24) 



Again the k is positive-definite. Note that for w — > oo the 
high-conductivity results ([T7| > are reproduced. In this case 
the k's are running with 1/r) while for w —¥ the k's 
are running with 1/r] 2 . This is a basic result: for low con- 
ductivity and for high conductivity the dependence of the 
k's on the magnetic Reynolds number Rm differs. For high 
conductivity (white noise) the k's are proportionate to Rm 
while for low conductivity (steep spectra) the factor Rm 2 
appears. Note that for <5-like spectral functions the numeri- 
cal coefficient for k is 0.2 while for k p this factor is about 
0.5. One can also find these values at the ordinate of Fig. [2] 

A basic difference exists for k and k p , too. While the n 
is positive-definite, the n p can change its sign. Generally it 
will be positive only for small w* but it should be negative 
for large w* . Already from these arguments one finds the 
main complication of the problem. The shape of the turbu- 
lence spectrum has a fundamental meaning for the results. 

To probe these results in detail a spectral function q(k) 



2L 



^ ~ 7T 1 + k 2 l 2 

is used. The integration yields 

1 St Rm v / w*(2 + 
15 
so that 



(1 + Vw*) 



^Rm St 

15 



for w* 



< 1 
> 4. 



(25) 



(26) 



(27) 
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Q_ 
< 
CL 
CL 

< 




Fig. 1 The k p vs. Rm after Eq. (128b . From top to bottom: 
St = 7, St = 4, and St = 1. All curves have a maximum. 
Note that all k p become negative for sufficiently large Rm; 
Pm = 1. 



Hence, for small w* (low conductivity) the k runs with 
St 5 Rm 15 while for high conductivity the relation is sim- 
ply St Rm. One finds again the differences between the two 
limits. In large-eddy simulations for the effective diffusivity 
the relation rj ~ ul c is used so that Rm ~ Pm ~ 1. As in 
the majority of the applications also the Strouhal number St 
is of the same order the coefficient (l27l i is a small number. 
If for direct numerical simulations the numerical value of 
Rm becomes large then there is no reason that 
smaller than unity. 



For k p there is another situation. One obtains 



1 RmSty / w*(3- 



15 



(1 



hence, 



iRm 15 St 5 

5 

-iRm St 



for 



< 9 
> 9. 



(28) 



(29) 



For w* > 9 the k p is negative so that the total pressure is 
always positive. For w* < 9, however, the k p becomes pos- 
itive. In this case for large Strouhal number the total mag- 
netic pressure 1 — k p becomes negative. Figure Q] demon- 
strates that K p exceeds unity for St > 4. For St > 4 one 
finds K p > 1 forRm = 14, i.e. w* — 3.5. We have to stress, 
however, that the SOCA approximation only holds if not 
both the quantities St and Rm simultaneously exceed unity. 
For turbulences in liquid metals in the MHD laboratory 
Rm ~ 1 is a typical value. Kernel et al. (2012) report an in- 
crease of k p with Rm 2 (their Fig. 9, bottom). The negative 
branch of d29l does not exist in the simulations (Pm = 0.5). 



3.2 Pm < 1 

The situation is more clear for small magnetic Prandtl num- 
bers, which exist, e.g., in stellar interiors, protoplanetary 
disks and also in the MHD laboratory. It is possible to con- 
sider the limit v — » in the Eqs. (fl4l but only for turbulence 



spectra with finite correlation time. Stationary patterns with 
E oc 5(ui) are excluded. In the limit of very small Pm the 
Eqs. (TBI) reduce to 



15t7 

c 

47T 

15r) 



E(k,Q) dk- 



E(k,0) dk 



1 

15 



E(k,u)k 2 



9 

15 



; 2 + r] 2 k 4 
E(k,uj)k 2 

LU 2 + T] 2 k 4 



dkduj, 



dkduj. (30) 



The spectrum (f2Tb leads to 

oc 

2-qk 2 



157710 J rjk 2 




1 



- w 
w 



8-qk 2 
15rjw J rjk 2 



q(k) dk, 



q(k) dk. 



(31) 



Again k is positive-definite. If the wave number spectrum 
has only a single value then 



Rm 2 2 



■ w 



15 



1 



_ 1 Rm 2 8 - w* 
Kp ~ 15 w* 1 + w*- K } 

The limit w* — !• is here not allowed. Again the k p is pos- 
itive (negative) for small (large) w* . Formally, the Strouhal 
number St does not appear. Replacing the w* by Rm/St in 
remains both limits the |k p | runs linearly with St Rm, i.e. with 1/rj. 

The k also runs with 1 jr\ in the high-conductivity limit, 



i.e. 



1 

15 



-St Rm. 



(33) 



while for low conductivity the value is n ~ (2/15)St Rm. 
For large w*, i.e. for Rm > St, and Pm <^ 1 there is prac- 
tically no influence of the numerical value of the magnetic 
Prandtl number (see Eq. |2~7| >. Below we shall also demon- 
strate by numerical solutions of the integrals that Eq. d33l 
forms the main result of the present analysis. Whether the 
K-coefficient may become larger than unity only depends on 
the numerical values of St and Rm. For large-eddy simula- 
tions with St = Rm = Pm = 1 the k is basically only of 
order 0.1. 

With the spectral function d25l l the results are very sim- 
ilar, i.e. 



1 



Rm 2 8 



15 



(34) 



1 + Vt»* 

This expression only exceed unity for St 3> 1. For small St 
the sum 1 — k p is thus always positive independent of the 
actual value of Rm contribution. 



4 Strong fields 

So far only the influence of weak magnetic fields has been 
considered. The influence of strong magnetic fields is also 
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10 



< 

< 




0.00 
-0.01 







10 



Fig. 2 



The /t/Rm 2 (top) and K p /Rm 2 (middle, bottom) 



vs. w* for the one-mode model ( l35l l. The curves in the plots 
(from fop fo bottom) are for S = 0.01, S = 1, S = 3, 
and S = 10. At the left vertical axis the values are valid 
for the delta function spectra (low-conductivity limit). The 
quantities vanish as 1/w* for w* — > oo (high-conductivity 
limit, right vertical axis) leading to the result ( 1331 ). Bottom: 
details for K p /Rm 2 ; Pm = 1. 



important to know. The rather complex results of the SOCA 
theory with arbitrary magnetic field amplitudes and with 
free values of both diffusivities are given in the Appendix. 
These expressions can be discussed by applying the single- 
scale wave number spectrum 



q(k) = 2u 2 5(k - l' 1 ) 



(35) 



and the frequency spectrum OTh . Such an approximation 
allows to solve the Eqs. (A1)...(A4) numerically including 




10 



w 




w 



Fig. 3 The same as in Fig.|2]for Pm = 0.1. From top to 
bottom: S = 0.01, S = 1, and S = 10. The quantities vanish 
as 1/w* forw* — > oo. 



the frequency integration so that the turbulence quantities 
ft/Rm 2 and K p /Rm 2 only depend on the Lundquist number 

S = (36) 



of the magnetic field, the frequency w* and the magnetic 
Prandtl number Pm. 

In the weak-field limit, S< 1, one finds the overall re- 
sult that k/Kto 2 runs as l/15w* (Figs. [2] and [3j top) so 
that again the general result (l33l is reproduced. For very 
small w*, i.e. for delta function frequency spectra (or, what 
is the same, for very long correlation times), the k's run with 
1 /Rm 2 - as already shown above. 

When the field is not weak, the stress parameters rapidly 
decrease with S. Figures [2] and [3] also demonstrate that the 
magnetic quenching can be written as 

K 



-1 + eS- (3?) 
(see Fig. |4|, in confirmation to Brandenburg et al. (2010) 
who found the magnetic quenching in terms of 1/B 2 . From 
the Figures one finds that e < 1 for Pm < 1. For large Pm 
the e is even smaller. The magnetic quenching of the k- 
parameter is thus stronger for small magnetic Prandtl num- 
ber than for large Pm. While a magnetic field with S = 1 
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reduces the n remarkably if Pm < 1 in the opposite case 
Pm > 1 the k is almost uninfluenced by S = 1. Figure 
^demonstrates the inverse dependence of the e on the mag- 
netic Prandtl number. One finds e ~ 0.75/Pm. The quench- 
ing expression, therefore, turns for Pm ^ 1 into 



1 + 0.75 Ha 2 



(38) 



with the Hartmann number Ha = S/ vPm instead of the 
Lundquist number S. For the magnetic quenching it is thus 
not important which of the diffusivities is large and which is 
small. The quenching is very strong if one of them is small 
(see Roberts & Soward 1975). For the high-conductivity 
limit (rj —> 0) or for inviscid fluids (v — > 0) the Hartmann 
number Ha takes very large values so that even very weak 
fields strongly suppress the K-effect. 
Note that 

S = Rm— , (39) 



B, 



cq 



with B cq = yl nop (u 2 ) as the equilibrium field strength. 
The magnetic quenching of the K-term thus grows with 
Rm 2 (Brandenburg & Subramanian 2005) so that for grow- 
ing Rm the k becomes smaller and smaller: 

1 St S 2 

K — (40) 

15e Rm B 2 

It becomes thus clear that in the high-conductivity limit even 
for rather small fields the K-term in Eq. (0 takes very small 
values which do not play an important role in the mean-field 
magnetohy drodynamic s . 

On the other hand, for Rm = 1 the well-known standard 
expression 

K 



1 



. B 2 



(41) 



for magnetic quenching appears with e of order unity only 
slightly differing for small and large w* . 

Because of St = Rm = 1 in this case the k's always 
remain smaller than unity in accordance to ( T33T >. Hence in 
both the possible concepts, i.e. the use of the microscopic 
diffusivities and the use of the large-eddy simulations with 
subgrid diffusivities, the values of the turbulence-induced 
Maxwell tensor coefficients remain small. 

The numerical simulations by Kernel et al. (2012) in- 
deed yield a magnetic quenching of the pressure term in 
terms of Rm 2 but only for Rm < 10. 

5 Catastrophic quenching? 

We have computed the stress tensor which is formed by 
large-scale background fields, by the Reynolds stress of a 
turbulence field under the influence of the field and the tur- 
bulent Maxwell stress of the field fluctuations. All contribu- 
tions can be summarized in form of the classical Maxwell 
stress tensor but with turbulence-modified coefficients (see 
Eq. [7}. The modified pressure term is now 1 — k p while 




10 



w 



Fig. 4 The verification of the relation ( 1381 for the func- 
tions w*k marked by their Lundquist numbers S. The re- 
sulting value for e is about 0.75; Pm = 1. 




0.1 



1 .0 
Pm 



10.0 



Fig. 5 The (weak) dependence of the quantity Pm • e on 
the magnetic Prandtl number Pm. 



the modified magnetic tension term is written as 1 — k. The 
quantities k and k p have been computed within the quasi- 
linear approximation (SOCA) which can be used if the min- 
imum of both the numbers St and Rm is (much) smaller 
than unity. As almost all turbulences fulfill the condition 
St ~ 1, the validity of SOCA requires Rm = ul c jj] < 1. 
Under this restriction the resulting k's are always smaller 
than unity. For all magnetic Prandtl numbers Pm we found 
k as positive so that the non-pressure force term {BV)B 
is reduced under the influence of turbulence. This is in par- 
ticular true for the coefficients of the angular momentum 
transport terms B^ Br and B^ B z which, therefore, be- 
come more and more ineffective in turbulent fluids. 

The sign of k p strongly depends on the magnetic Prandtl 
number. It proves to be negative-definite for large Pm. For 
smaller Pm the sign of k p depends on the shape of the fre- 
quency spectrum of the turbulence. For steep profiles, i.e. 
very long correlation times, the k p becomes positive while 
for flat frequency-spectra of the turbulence which are as flat 
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as the spectrum of white noise (very short correlation times) 
the k p for Rm < 1 becomes negative. 

One could believe that relations valid for small Rm like 
k p cx St • Rm can be also used for Rm > 1 so that finally the 
effective magnetic pressure becomes negative. This, how- 
ever, is not true. The k p changes its sign for Rm 3> St and 
becomes negative. Hence, the total magnetic pressure re- 
sults as mostly positive. The only exception exists for suffi- 
ciently large St and sufficiently small Rm (see Fig.Q]). 

More dramatic is the situation with the magnetic ten- 
sion and its coefficient 1 — k which is also the coefficient of 
the vector J x Bin the generalized Lorentz force in turbu- 
lent media. This coefficient is positive for small k, i.e. for 
sufficiently small Rm if St = 1. It is positive and smaller 
than unity for the large-eddy simulations ('mixing-length 
model') considered at the end of Sect. l3.2l with Rm = St = 
Pm = 1 (see Fig. [2J. 

The question, however, whether the k can exceed unity 
(so that 1 — k becomes negative) cannot finally be answered 
within the quasilinear approximation. It is k ~ 0.1 St • Rm 
where one of the factors St and Rm must be smaller than 
unity but the product St • Rm is formally not restricted by 
the SOCA. It is thus a clear and surprising result also in 
the frame of SOCA that the angular momentum transport 
by large-scale magnetic fields can strongly be suppressed 
under the influence of turbulence. The possible existence of 
an instability resulting from n > 1 has been confirmed by 
the numerical simulations by Brandenburg et al. (201 1). 

The formal background of this phenomenon is that 
the integrals defining k and n p do not exist in the high- 
conductivity limit or, what is the same, in the ideal MHD. 
The same is true for the much simpler magnetic-suppression 
problem of the eddy diffusivity. We take the expression 



V k 2 E 

,2 _|_ ^2£,4 



^ 6 7] 2 k 4 



B z 



dfcdu; 



(42) 



5 (uj 2 + rj 2 k 4 ) 2 Hop, 

(Kitchatinov et al. 1994) for the SOCA expression of the 
eddy diffusivity under the presence of a uniform magnetic 
background field (Pm = 1). The expression is part of a se- 
ries expansion which converges if the second term is smaller 
than the first term. The second term of the RHS of this ex- 
pression has two important properties: i) it is positive for all 
spectral functions E with dE / du < so that the i]t is al- 
ways reduced by the magnetic fields, and ii) it does not exist 
for the limit r\ —> 0. In other words, for rather small T) the in- 
tegral becomes large so that the magnetic quenching would 
be extremely effective for large Rm. This is why such a 
series expansion only holds for very weak fields. This phe- 
nomenon has been called a 'catastrophic' quenching (see 
Blackman & Field 2000; Blackman & Brandenburg 2002). 
It exists within the SOCA theory for the eddy diffusivity 
and also for the eddy viscosity. One finds from Eq. d42l that 
the mentioned diffusivities are magnetically quenched like 



1 — S 2 for small S and like S~ 3 for large S. Of course, by 
this procedure the ryx cannot become negative. We know, 
on the other hand, that the magnetic quenching of the eddy 
diffusivity in sunspots reduces its value (only) from 5x 10 12 
cm 2 s _1 to about 10 11 cm 2 s _1 what - together with the 
time decay law of the sunspots - can be understood with 
quenching expressions like (l42b for Rm = 1 (Riidiger & 
Kitchatinov 2000). It is thus suggested to work with the sim- 
ple relations Rm = 1 and S ~ B /B cq in applications with 
turbulent convection. 

Similarly, also the n increases for vanishing rj. There 
is, however, no nonmagnetic term against which the mag- 
netic influence can be neglected as it must be compared 
with the large-scale Lorentz force J x B which is also of 
the second order in B. The only possibility to keep the tur- 
bulence contribution small for large Rm is to put St <C 1. 
However, if the magnetic field is super-equipartitioned then 
the k is magnetically quenched which introduces a new 
factor Rm -2 . Then the magnetic-induced K-effect finally 
runs with 1 /Rm so that it vanishes in the high-conductivity 
limit. In summary, for large Rm and for very weak mag- 
netic field the k can exceed unity (so that the stress tensor 
reverses sign) but this phenomenon disappears already for 
rather weak fields. 
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A SOCA expressions of the k's 

The expressions for the mean-field Lorentz force parameters k and 
ftp of Eq. 10 provided by the quasilinear theory for arbitrary mag- 
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netic amplitudes can be written as 





OO CO 



g^*p(*,*,*)dfafc,. (Al) 2 ^ COS ^ 4/34 V ^^sinf 

to -\- rj k A R \ 

-(/3 2 + l + 2cos0)— . (A4) 
The kernel functions AT and Jf p depend on the magnetic field and M 2 / 

the variables fc and u via The first parts in these expressions represent the contribution of the 

_ kV Reynolds stress while the following lines represent the small-scale 

~ (lu 2 +rj 2 k 4 y/ 4 (uj 2 + v 2 k 4 ) 1/4 ' Maxwell stress. 

B 2 - 2/3 sin | + T 



LN = log 



/3 2 + 2/3 sin § + 1 



/ B — sin § \ / /3 + sin J „ , 
Ai? = arctan + arctan —2- . (A2) 



Here, cos0 = (r,vk A - lu 2 ) /^/(uj 2 + ri 2 k 4 )(v 2 k 4 + lu 2 ). 
The kernels read 

U 2 + ^ 2fc V V J 4/3sinf + 

(^_ 3) ^fl \ 1 / 2 _ 

; 2/3cos|y 8/3 4 ^ VP ^4/3sin| 

2/3 cos | 



A R 

- (/3 2 +3 + 6cos</>) — -I (A3) 



and 
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